UH&AMtSi 


ELECTS 

OCT  1  8 1984 


AMTE  tTfiWteftan) 
9ue«n's  ftoa*  TEf&MfiTQN 

MkMtuix  twuow 


UNLIMITED 

AMTE(N)TM84058 


J.H.  JAMES 


FLUID  LAYER  BETWEEN  INFINITE  ELASTIC  PLATES 
II.  DISTRIBUTION  OF  POWER  FLOW 


BY 


|  Accession  For 

NTIS  GRAM 
DTIC  TAB 
Unannounced 
Justificati< 

jSSK'flS! 

■ 

Bv 

Distribution/ 

Availability  Codes 

Dist 

AlL 

Avail 

Spec 

Formulae  axe  given  which,  together  with  formulae  contained  in 
Part  I,  enable  numerical  evaluation  of  the  separate  power  flows  in  the 
plates  and  fluid  due  to  line-force  and  line-source  excitation.  Plots 
of  power  flow  show  significant  interchange  of  energy  between  plates  and 
fluid.  They  are  complimentary  to  the  acoustic  intensity  vector  plots 
of  Pa rt  I. 


14  pages 
5  figures 


AMTE( Teddington ) 

Queen's  Road 

TEDDINGTON  Middlesex  TW11  OI/I  March  1984 

© 


Copyright 

Controller  HMSO  London 


UNLIMITED 


A  previous  report  [l]  has  given  the  formulas  necessary  for  eval- 
uating  numerically  fha  pxaaauxa  and  displacements,  dua  to  line-fore*  and 
lina-aourca  excitation,  of  a  ayataa  which  comprises  of  two  identical 
ataal  plate*  separated  by  a  layer  of  fluid.  Plot*  of  intensity  vectors, 
in  a  20c*  layer  of  water  bounded  by  lea  steal  plates,  showed  significant 
interchange  of  energy  between  the  fluid  and  the  driven  plate  When  the 
excitation  was  a  line-force,  but  not  When  the  excitation  was  a  source. 

One  of  the  recoanwnd&t ions  of  the  report  was  that  the  distribution  of  the 
power  flow  between  plates  and  fluid  should  be  evaluated  as  a  function  of 
the  horizontal  distance  from  the  excitation,  this  is  the  main  aim  of 
the  present  report. 

The  motivation  behind  this  suit*  of  reports  is  the  need  to  under¬ 
stand  the  physics  of  the  interchange  of  energy  between  a  pipe's  wall  and 
its  contained  fluid.  James  [2]  gives  intensity  vector  plots  Which  shear 
this  effect  in  a  steel  pipe  containing  water,  and  Puller  [3]  gives  some 
formulae  Which  enable  the  distribution  of  power  flew  between  the  pipe's 
wall  and  the  fluid  to  be  calculated,  at  large  distances  fro*  the  source. 

The  distribution  of  power  varies  with  distance  due  to  the  non-orthogonality 
of  the  system  modes  over  the  fluid's  cross-section.  At  large  distances 
fro*  the  source,  it  is  the  interaction  of  predominantly  'fluid'  modes  with 
predominantly  'elastic*  sodas  which  is  responsible  for  the  interchange  of 
energy.  While  close  to  the  source,  interaction  among  evanescent  wave-fields 
may  cause  intense  energy  circulation  between  the  fluid  and  wall. 


2.  MATHEMATICS  OP  PROBLEM 


The  pressure  and  displacements  of  the  syste*  shown  in  Figure  1  are 
represented  by  the  Fourier  integrals 
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in  which  Wj(x)  and  W2(x)  are  the  displacements  of  the  upper  and  lower  plates 
respectively/  p(x,s)  is  the  acoustic  pressure  in  the  fluid  layer/  Mg(x,s) 
is  the  horizontal  acoustic  particle  displacement .  The  time-harmonic  factor, 
*iqp(-iwt),  is  oaiitted  throughout.  Formulae  for  the  transforms  Wj{a)  etc. 
are  given  elsetdiere  [1]. 


With  reference  to  Figure  18,  which  shows  a  sign  convention  for  the 
positive  directions  of  Shear  forces  and  moments,  the  power  per  unit* length 
flawing  in  the  positive  x-dlrectlon  of  each  plate  is  the  sum  of  the  two 


(2) 


PM(x)  -  ( 1/2 )Beal[M( x ) . -dw  (x)/dx] 
*s(x)  -  (1/2 )Real[-S(x)W*(x)] 


In  each  platn.  p^x)  and  Fs(x)  an  the  powers  transmitted  by  the  bending 
moments  and  shearing  forces,  respectively;  W(x)  is  the  normal  velocity  of 
the  plate,  -1 wN(x)i  and  the  symbol  *  denotes  complex  conjugate.  The 
bending  moment  and  the  Shearing  stress  per  unit  length  are  simply 

^  M(X)  -  D.d2t»(x)/dx2  ( 3 

s(x)  -  -o.d3w(x)/dx3 

where  0  Is  the  plate's  flexural  rigidity.  The  following  Fourier  integrals 
are  obtained  from  equations  ( 1 )  and  ( 3 ) 
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The  power  flow  per  unit  length  in  the  fluid,  in  the  positive  x- 
direction,  is  the  integral  of  the  horizontal  component  of  the  acoustic 
intensity  over  the  fluid's  cross-section.  It  is 

H 

Pp<x)  -  ( 1/2 )  |  Real[p< x, z )W*( x, z ) ]dz  (5) 

0 

In  the  case  of  line-force  excitation,  F26(x),  a  check  on  the 
numerical  results  is  possible  by  use  of  the  formula 

PM1<X)^S1(X)+PM2(X>+PS2<X)+PF<X>  “  <V4>Real[F2**(0)J  (6) 

which  is  only  valid  close  to  x-0,  because  dissipation  is  included  in  the 
system.  Equation  (6)  reflects  the  fact  that  the  input  power  of  the  line- 
force  excited  plate  flows  equally  in  the  positive  and  negative  directions. 


3.  WHMOUCM.  RESULTS 

The  Fourier  integrals  of  equation  (4),  with  upper  limits  of  integ¬ 
ration  set  to  approximately  twice  the  highest  fres  wavenumber  at  the 
selected  frequency,  were  evaluated  by  a  simple  adaptive  Oausslan  quadrature 
scheme.  Because  the  integrals  must  be  evaluated  in  the  principal  value 
sense,  it  is  necessary  to  introduce  damping  into  the  system  via  a  complex 
Young's  modulus,  E(l-ir)g),  and  a  complex  sound  velocity,  c(l-it>f).  The 


Integral  In  equation  (5)  was  evaluated  by  Simpson's  rule,  from  the  21 
values  of  p(x,z)  and  Wx<  x,  z )  that  were  obtained  at  each  of  60  horizontal 
x- st art: ions.  At  x-0,  the  power  transmitted  by  the  shearing  force,  S(x), 
cannot  be  evaluated  numerically,  but  evaluation  proceeds  with  increasing 
accuracy  as  x  increases;  the  power  transmitted  by  the  bending  moment  is 
zero  by  virtue  of  dW/dx  vanishing;  and,  the  acoustic  power  flow  must  be 
zero  because  Wx  vanishes,  subject  to  the  accuracy  of  the  computations, 
the  consistency  test  of  equation  (6)  was  satisfied  by  the  numerical  results 
shown  herein,  the  following  constants  in  SI  units  were  used  to  obtain 
Figures  2-5 t 

E-19.5X1010  ff-0.29  p^-7700.0  h-0.01 

p-1000.0  0-1 500.0  H— O .20  Z0»H/3 

n  <-0.02  n -“O-ooi 
s  f 

The  dispersion  plots  of  the  dissipation-free  system  are  shown  in 
Figure  2.  The  near  identical  branches  labelled  1  and  2  are  symmetric  and 
antisymmetric  'plate*  waves  whose  energy  is  mostly  confined  to  the  plates, 
their  group  velocity  being  approximately  twice  their  phase  velocity.  The 
branches  labelled  3-5  are  predominantly  ’fluid'  waves  in  which  the  phase 
and  group  velocites  at  cut-on  are  infinity  and  zero,  respectively. 

In  Figures  3-5  the  percentages  of  the  total  power  flow,  in  the 
positive  x-direction,  contained  in  the  plates  and  fluid  is  given  as  a 
function  of  distance  from  the  excitation.  Power  is  not  conserved  because 
of  finite  values  of  the  loss-factors,  but  the  attenuation  is  small  - 
in  fact,  not  more  than  ldB  in  the  range  x-0. 4  to  0.9m.  The  power  in  the 
top  plate  is  very  small  compared  with  the  power  in  the  bottom  plate,  as 
can  be  inferred  from  the  intensity  vector  plots  [1]. 

The  power  flow  due  to  a  line-force  excitation  is  shown  in  Figures 
3  and  4.  At  lkHz  the  individual  power  flows  settle  down  to  constant  values 
due  to  the  absence  of  'fluid'  waves,  the  power  in  the  fluid  being  less  than 
10%  of  the  total.  At  4,  7  and  10kHz  the  power  levels  away  from  the  source 
oscillate  considerably  due  to  interaction  between  'plate'  waves  (1,2)  Whose 
wavelengths  are  almost  equal,  and  'fluid'  waves  (3,4,5)  of  differing  wave¬ 
lengths  .  At  4kHz  the  power  levels  oscillate  almost  sinusoidally  because 
there  is  only  a  single  'fluid*  wave,  the  level  in  the  fluid  ranging  from  a 
minimum  of  10%  to  a  maximum  of  30%;  at  7kHz  the  presence  of  two  'fluid' 
waves  causes  complex  beating  in  which  the  level  varies  from  15%  to  45%; 
and,  at  lOkhz,  where  three  'fluid'  waves  have  cut-on,  the  power  in  the 
fluid  ranges  from  15%  to  50%.  In  the  aforementioned  plots,  both  the  fluid 
and  plate  power  flows  are  positive  everywhere.  However,  the  high-amplitude 
evanescent  wave-fields  near  to  cut-on  frequencies  sometimes  superpose  local 
circulatory  energy  flows  between  fluid  and  plate  Which  result  in  net  local 
power  flows  in  the  fluid  Which  may  be  negative.  Some  plots  obtained  at 
2  and  5kHz,  but  not  included  here,  show  this  effect. 

Intense  circulating  energy  flows  between  plates  and  fluid  have  been 
illustrated  elsewhere  [1]  in  a  plot  of  Intensity  vectors  due  to  a  line- 
source  excitation  at  1kHz.  The  power  flow  is  shown  in  the  top  of  Figure  5. 
The  local  power  flow  in  the  fluid  ranges  from  flow  in  the  negative  direction 
of  190%  to  flow  in  the  positive  direction  of  66%,  While  the  power  in  the 
plate  ranges  from  280%  to  34%  of  the  net  power  flow.  At  distances  greater 


V 


than  0.8a  from  th*  souse*,  th*  levels  hav*  settled  down  to  th*  mm  constant 
jpsrnsntsgss  as  wszs  obtained  fxoa  line-fore*  excitation.  At  £c*qu*ncl*s  of 
4,  7  and  lOkBx,  th*  power  in  th*  plats*  dus  to  line-source  switatlon  is 
easy  small;  hones,  only  th*  plot  at  «cBs  is  Shown,  in  th*  lower  half  of 
Figure  5. 


4. 


Pornula*  hav*  boon  given  which  snabls  numerical  svaluatlon  of  the 
separate  power  flows  in  the  plats*  and  fluid.  Plots  of  pow*r  flow  have 
b*sn  presented  for  the  ease  of  a  20ea  layer  of  water  bounded  by  lot  steel 
plates,  they  add  quantitative  information  to  the  qualitative  information 
contained  in  the  intensity  vector  plots  given  elsewhere  [1].  When  the 
excitation  is  a  line-foxes  the  plots  (a)  Shew  that  most  of  the  power  flow 
resides  in  the  plates  ;  (b)  help  to  confirm  that  the  energy  interchange 
is  caused  by  an  interaction  between  'plate'  and  'fluid*  waves >  (c)  Shew 
that  inergy  interchange  nay  be  significant  -  up  to  35%  of  the  total  at 
lOkHs.  When  the  excitation  is  a  line-source,  the  plots  (a)  Show  negative 
power  flow  close  to  the  source,  due  to  intense  circulating  energy  flow 
caused  by  interaction  among  evanescent  waves i  (b)  show  negligible  power 
in  the  plates  tftien  'fluid'  waves  have  cut-on. 

Future  reports  will  extend  the  rang*  of  material  and  geometric  con¬ 
stants,  and,  in  particular  will  (a)  Show  plots  at  frequencies  close  to 
the  cut-on  frequencies  of  'fluid'  waves  and  (b)  investigate  energy  inter¬ 
change  when  the  distinction  between  fluid  and  plat*  waves  is  less  clear. 
■  Also  to  be  investigated  is  the  effect  of  line-constraints. 
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